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ABSTRACT 

We compare models for Type la supernova (SN la) light curves and spectra with an extensive 
set of observations. The model s come from a recent survey of 4 4 two-dimensional delayed- 
detonation models computed by lKasen. Ropke & Woosleyl (l2QQ9l) . each viewed from multiple 
directions. The data include optical light curves of 251 SNe la, some of which have near- 
infrared observations, and 2231 low-dispersion spectra from the Center for Astrophysics, plus 
data from the literature. These allow us to compare a wide range of SN la models with observa- 
tions for a wide range of luminosities and decline rates. The analysis uses standard techniques 
employed by observers, including MLCS2k2, SALT2, and SNooPy for light-curve analysis, 
and the Supernova Identification (SNID) code of Blondin & Tonry for spectroscopic compar- 
isons to assess how well the models match the data. The ability to use the tools developed for 
observational data directly on the models marks a significant step forward in the realism of 
the models. We show that the models that match observed spectra best lie systematically on 
the observed width-luminosity relation. Conversely, we reject six models with highly asym- 
metric ignition conditions and a large amount (> 1 Mq) of synthesized ^^Ni that yield poor 
matches to observed SN la spectra. More subtle features of the comparison include the gen- 
eral difficulty of the models to match the [/-band flux at early times, caused by a hot ionized 
ejecta that affect the subsequent redistribution of flux at longer wavelengths. The models have 
systematically higher velocities than the observed spectra at maximum light, as inferred from 
the Si II 16355 line. We examine ways in which the asymptotic kinetic energy of the explosion 
affects both the predicted velocity and velocity gradient in the Si II and Ca II lines. Models 
with an asymmetric distribution of ^^Ni are found to result in a larger variation of photometric 
and spectroscopic properties with viewing angle, regardless of the initial ignition setup. We 
discuss more generally whether highly anisotropic ignition conditions are ruled out by obser- 
vations, and how detailed comparisons between models and observations involving both light 
curves and spectra can lead to a better understanding of SN la explosion mechanisms. 
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1 INTRODUCTION and precise cosmo logical distances t o establish the acceler ation of 

, . , cosmic expansion (iRiess et aPI 19981 : IPerlmutter et alJll999h . 

Type la supernovae (SNe la) play a major role m many astrophys- 

ical pheno mena. They produce a large fraction of iron in the uni- Respite their astrophysical importance, however, they remain 

verse (e,g.,'Truran & Camerodllizll), heat the interstellar medium enigmatic objects. There is a general consensus that they result 

(e.g.,|ciotti et al. 1991), and form an endpoint of binary star evolu- f™™ *e thermonuclear disruption of a carbon-oxygen white dwarf 

tion fe.g.. llben&Tutuk^[T984h . SNe la provide the most reliable ^^D) star Moyle & Fowle Jll960D approaching the Chandrasekhar 

mass (Mch ~ 1.4 Mq), either through accretion from a non- 
degenerate binary companion (the "single-degenerate" scenario), 
or thr ough merger with another WD (the "do uble-degenerate" sce- 
* E-mail: blondin@cppm.in2p3.fr nario; llben & Tutukovlll984l : IWebbink|[T984h . Which of these two 
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possibilities constitutes the d ominant (or so le) progenitor channel 
for SNe la is still debated fsee lHowelfcOld for a recent review). 

The explosion mechanism it self is also largely unknown 
(see iHillebrandt & Niem ever 200(j for a review). In the preferred 
"delayed-detonation" model (Khokhlov 1991), the burning starts 
as a turbulent subsonic deflagration near the WD center and 
transitions to a supersonic detonation near its surface. The de- 
flagration pre-expands the WD so that the subsequent detona- 
tion does not burn the entire star to nuclear statistical equilib- 
rium (NSE) material (including ^^Ni to power the light curve), 
but instead synthesizes appropriate fractions of high-velocity (^ 
10000 km s~^) intermediate-mass elements (IME; such as Mg, 
Si, S, Ca, etc.) needed to reproduce the observed spectra. The 
first simulations were carried out in ID, but recent studies show 
that multi-dimensional simula tions are needed to capture hy- 
drodynamical instabilities (e.g. jGamezo. Khokhlov & Oranll2005l : 
iRopke & Niemevej|2007h and to provide a physi cal basis for the 
transition from deflagration to d etonation, (e.g., IWooslevI 120071 : 
lRopkell2007l :IWooslev et al." 2009). 

The empirical relation between the peak luminosity and the 
width of the light-curve (the so-called w idth-luminosity relation, 
or WLR: |Pskovskiilll97l |Phillipslll99 3'). instrumental to the use 
of SNe la as distance indicators, can be physically interpreted in 
terms of (1) varying opacity with the amount of synthesized ^^Ni 
(iHoeflich & Khokhlov! 1 19961). (2) varying mass of the progenitor 
WD (|Pinto&Eastmanll2000ah . or (3) the lll^ll (i.e. doubly- to 
singly-ionized) recombination timescale of iron-group elements in 
the SN atmosphere ( Kasen & Woosley 2007). More detailed calcu- 
lations are also needed to exp lain the observed sc atter in the WLR. 

In a recent paper, .Kasen, Ropke & WooslevI (.2009.) (hereafter 
KRW09) conducted a 2D survey of delayed-detonation models, in 
which they varied the radial/angular distribution and number of ig- 
nition points in a Chandrasekhar-mass C-i-0 WD star, as well as the 
criterion for deflagration-to-detonation transition. By synthesizing 
light curves and spectra for different viewing angles, they were able 
to broadly reproduce both the observed width-luminosity relation 
(for all but subluminous SNe la) and its scatter, illustrating the im- 
portance of multi-dimensional computations to reproduce observed 
trends in SN la properties. Furthermore, they showed that variations 
in the metallicity of the progenitor WD affect both the slope and 
normalization of the WLR, and that ignoring these effects could 
potentially lead to systematic overestimates of '-^ 2% on distance 
determinations to SNe la. 

To study such subtle effects, one must ensure the models re- 
produce all observed trends in some detail, and not just the width- 
luminosity relation. This is precisely what we set out to do in 
the present study, where we conduct an extensive analysis of the 
delayed-detonation models of KRW09 through a detailed and di- 
rect comparison with observations of SNe la. We present the mod- 
els and data in § [3 and our methods for quantitatively evaluating 
each model in § [3] We then proceed to a detailed comparison of 
their photometric (§ |4]) and spectroscopic (§ [5]) properties with ac- 
tual data. We discuss whether explosion models with asymmetric 
ignition conditions are ruled out by observations in § [6l and con- 
clude in §[71 



dimensional delayed-detonation explosions in a non-rotating spher- 
ical Chandrasekhar-mass WD star composed of equal mass frac- 
tions of carbon and oxygen. They varied both the radial/angular 
distributions and number of ignition points used to trigger the de- 
flagration (which control the level of asymmetry of the explosion), 
as well as the criterion for the transition to a detonation (termed 
"dc", and parameterized via the critical Karlovitz number, which 
controls the level of turbulence in the combustion). We refer the 
reader to KRW09 for more information on the initial conditions for 
each model, but note that these were chosen a priori, and not tuned 
to match the range of observed SN la properties. 

Table [T] gives the asymptotic kinetic energy (^kin), abun- 
dances, peak bolometric luminosities (Lboi,peak) and Ami5(Bjj 
decline-rate ranges for all the models of KRW09. Models with 
an isotropic (anisotropic) distribution of ignition points are la- 
beled DD2D_iso (DD2D_asym). One exception concerns the 
DD2D_asym_01 model series, which in fact has an isotropic dis- 
tribution of ignition points. The number that follows (01-08) 
corresponds to a particular ignition setup, and the criterion for 
deflagration-to-detonation transition (dc) is also included in the 
model name (see Tables 1 & 2 in KRW09). The models span a 
large range in mass of synthesized ^^Ni (0.29-1.10 Mq) and de- 
cline rate parameters Ami5{B) (0.74-1.57 mag). The models are 
unable to reproduce rapidly-de clining SNe la such as the sublu- 
minous SN 1991bg fsee .Taubenberger et alJuO OS for a discussion 
on this class of objects). Model DD2D_iso_01_dc4 was accidently 
omitted in the radiative transfer calculations, but this does not affect 
the results presented in this paper. We do not study the impact of 
metallicity variations in the progenitor WD star (between 1/3 to 3 x 
solar; see KRW09). We selected eight models for detailed analysis, 
and rejected six models based on their spectra (see § I3.31 l. These 
are marked as "subset" and "rejected" in Table[T] 

The models were evolved hydrodynamically until ^ 100 s 
past ignition, at which point the ejecta had reached a phase 
of homologous expansion (velocity proportional to radius). Fig- 
ure [T] shows the total mass density distribution for models 
DD2D_iso_06_dc2 (top left) and DD2D_asym_01_dc3 (bottom left). 
The white contours correspond to a ^^Ni mass fraction of 10~^. 
The ^^Ni distribution is roughly symmetric about (z = 0) in model 
DD2D_iso_06_dc2, while it is strongly skewed toward z > in 
model DD2D_asym_01_dc3, this despite the isotropic distribution 
of ignition points in both models. 

A crude nuclear network was used to determine the distribu- 
tions of radioactive ^^Ni, stable iron-group elements (IGE), and 
intermediate-mass elements (IME; such as Mg, Si, S, and Ca), but 
detailed abundances of all important chemical species were deter- 
mined based on a more elaborate nuclear network applied to a rep- 
resentati ye model (DD2D_iso -06-dc43) using 12800 tracer parti- 
cles f see lTravaglio et al.l2004h and interpolated accordingly for the 
other models. The nucleosynthetic yields are thus subject to some 
uncertainty which can impact the derived observables. 

The abundances were then remapped on a 64 x 128 cylin- 
drical grid with a velocity cutoff of ^ 25000 km s~^ for 
the radiative transfer calculations. These were done with the 
time-dependent Monte Ca rlo radiative transfer code SEDONA 
(IKasen. Thomas & Nugent 2006) using 10^ photon packets and 13 
chemical species (ionization stages I-V): C, O, Na, Mg, Si, S, Ar, 



2 MODELS AND DATA 

2.1 2D delayed-detonation models from KRW09 

iKasen. Ropke & WooslevI (120091) simulated 44 axisymmetric two- 



^ the difference in B-band magnitude between maximum light and 15 d 

after maximum; Ph illipsiil993 ) 

2 and not DD2DJso_06_dcl, as stated by KRW09. 
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Figure 1. Top left: Total mass density distribution at t = 100 s for model DD2D_iso_06_dc2. The ejecta are in homologous expansion, and 10^^ cm in 
(r, z) corresponds to 10^ km s"-*^ in (^'r, Vz) . The white contours correspond to a ^^Ni mass fraction of 10~^. The arrows denote the 30 different viewing 
angles. Top center: Bolometric light curves as a function of viewing angle. The inset shows the maximum-light region, with filled circles corresponding 
to the time of bolometric maximum. Top right: Optical spectra (3000-8000 A) as a function of viewing angle. Bottom row: Same as above but for model 
DD2D_asym_01_dc3. 



Ca, Ti, Cr, Fe, Co, and Ni. Local thermodynamic equilibrium (LTE) 
was assumed for the atomic level populations, although non-LTE 
effects in the radiation field were treated approximately through 
an equivalent two-level atom formalism. The photon packets were 
collected into 30 separate viewing angle bins (15° < < 165°), 
offering a 2D view of the explosion. The bins are equally spaced 
in cos 0, so that each bin subtends the same solid angle and has an 
equal probability of being observed. The arrows in the left panels 
of Fig. [T] denote the viewing angles considered by KRW09. The 



limited number of photon packets and the large number of view- 
ing angles results in moderate S/N in individual angle-dependent 
light curves and spectra (typically S/N^ 15-20 per 10 A in B and 
V at maximum light). The velocity cutoff of ^^ 25000 km s~^ also 
limits the formation of high- velocity absorption features, which are 
common in early-time SN la spectra (e.g. Mazzali et al. 2005). 

Spectra were computed between 2.5 d and 44.5 d after explo- 
sion (i.e. between — 17( — 15) d and -i-25(-i-27) d from B-band max- 
imum for a 19.5(17.5) d rise time) in 1 d steps. At later times, non- 
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LTE effects and non- therma l excitation by fast electrons become 
increasingly important (e.g., Kozma & Franssonlll992 ) and the ra- 
diative transfer calculations are less reliablo The mass fraction 
of (unburnt) carbon dominates over all other intermediate-mass el- 
ements at V > 1 6000-2 1000 km s~^, but we are unable to study 
the presence of corresponding spectral lines (e.g. C II 16580) in the 
early-time spectra due to limited signal in the Monte Carlo spectra. 
The total mass of (unburnt) carbon is < 0.04 Mq for all models. 

Figure [T] shows the variation of bolometric light curves {mid- 
dle panels) and optical spectra {right panels) with viewing angle for 
models DD2DJso_06_dc2 {top) and DD2D_asym_01_dc3 {bottom). 
In model DD2D_iso_06_dc2, the peak bolometric luminosity varies 
by < 10% with viewing angle, and the spectra only change sig- 
nificantly blueward of ^ 3500 A. In model DD2D_asym_01_dc3, 
however, the peak bolometric luminosity varies by ^ 70% with 
viewing angle, and the spectra are affected over a large fraction of 
the optical range (out to ^ 6000 A). 

These variations are not related to the total mass of ^^Ni 
synthesized during the explosion (^O.TOMq and ^0.64 Mq for 
DD2D_iso_06_dc2 and DD2D_asym_01_dc3, respectively); rather, 
they are a consequence of the distribution of ^^Ni in the ejecta. 
When the model is viewed from directions where ^^Ni extends to 
larger radii, the bolometric light curves peak earlier at a higher lu- 
minosity and are broader than for viewing angl es where ^^Ni is 
confined to the deeper ejecta regions (see also IPinto & EastmanI 
l2000al . their Fig. 4). The impact on the spectra is dramatic, ow- 
ing to the higher ejecta temperatures (resulting in a bluer SED) 
and higher ionization (affecting the relative shapes and strengths of 
spectral features). For both models, the variation is largest at ultra- 
violet wavelengths, and is related to differences in the abundance 
of iron-group ele ments in the outer layers of the ejecta (see, e.g., 
ISauer et al.ll2008h . 

The variation in peak bolometric luminosity with viewing an- 
gle in a given model (see Table [B can thus be used to gauge the 
level of asymmetry in the ^^Ni distribution, but it does not constrain 
the level of isotropy in the distribution of ignition points. Most 
DD2D_iso models display variations in Lboi,peak with viewing an- 
gle well below the 20% level, but others (DD2D_iso_04_dc3 and 
dc4) display factor-of-two variations. Likewise, some DD2D_asym 
models with the highest level of anisotropy in the distribution 
of ignition points (e.g. DD2D_asym_07_dc2 and dc3) display < 
10% variations in peak bolometric luminosity with viewing angle, 
whereas the DD2D_asym_01 model series (which has an isotropic 
distribution of ignition points; see above) displays the largest vari- 
ations (65-70%). This is not entirely surprising, since the ignition 
points are randomly distributed within a given setup, and individual 
ignition points that happen to placed farther out than the bulk of the 
ignition sparks can easily dominate the flame morphology (this is 
especially true in 2D, which tend s to favour large Rayleigh-Taylor- 
structures; see lRopke et alj|2006h . 



2.2 Data 

We largely rely on our database of light curves and spectra from 
the CfA Supernova Program to compare the models w ith observa- 
tions. These include opt i cal light curves of 25 1 SNe la JRiess et alj 
Il999bl : IJha et all l2006l : iHicken et alJ 120091) . with near-infrared 
(NIR) photometry for the brighter objects dWood-Vasev et alj 
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Figure 2. Distribution of the Amis (B) decline rate parameter in our SN la 
sample with a maximum-light spectrum. SNe la with spectra resembling 
those of the luminous SN 1991T and the faint SN 1991bg are highlighted. 
Names of peculiar SNe la are given in their respective Ami5(B) bin. 
The horizontal error bar shows the Amis (5) range in the 2D models of 
KRW09. 



|2008% and 2231 optical spectra o btained for th e most part us- 
ing the FAST spectrograph (Fabricant et al.lll998h mounted on the 
1.5 m Tillinghast telescope at the Fred Lawrence Whipple Observa- 
tory (FLWO). The spectra have a typical FWHM resolution of 6- 
7 A with a rest-frame wavelength range 3500-7500 A. Already 577 
SN la spectra have been published in several papers (the largest 
collection of 432 sp ectra of 32 SNe la has been published by 
iMatheson et al.ll2008h and are publicly available via the CfA Su- 
pernova ArchivO Additional spectra will be published in a forth- 
coming paper. 

We complemented this data set with published data from 
the literature, and reference our sources where appropriate. Our 
study of rise times (§ 14.11 ) makes use of measurem ents on SN la 
light curves from the SDSS-II Supernova Survey bv iHavden et al.l 
(2010), as well a s pseudo-bolometric light curves published by 
IStritzingeJ(l2005h . 

The observed SN la sample used in this paper offers a fair 
representation of the true SN la population in the local Universe. 
Figure |2] shows the Ami5{B) distribution of SNe la for which 
we have a spectrum within three days from maximum light. This 
sample spans a large range in Ami5{B) (^0.7-2.1 mag) and in- 
cludes fractions of luminous 1991T-like (^ 9%), faint 1991bg-like 
{^ 7%), and pecuhar 2002cx-like (-> 3%) SNe la comparable to 
those found by Li et aL| (1201 Ih for a volume-limited sample (9%, 
15% and 5%, respectively). The fraction of 1 991bg-like SNe la is 
a factor of two smaller than that reported by iLi et al.l (1201 ih and 
reflects a possible magnitude bias (the fraction drops from 15% 
to 3% for an ideal magnitude-limited sample; Li et al. 2011, their 
Fig. 11). In any case, the models of KRW09 do not extend beyond 
Amis (5) ~ 1.6 mag and none of their spectra present the promi- 
nent Ti II absorption band around 4000-4500 A characteristic of 
1991bg-like SNe la. The fraction of 2002cx- like SNe la is also 
smaller than for the volume-limited sample of iLi et all (l201lh . but 
the reported rate has a large associated error. 



^ We do note, however, that KRW09 present light curves for model 
DD2D Jso_06_dc2 until +60 d past S-and maximum in their Fig. 2. 



^ http://www.cfa.harvard.edu/supernova/SNarchive.html 
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3 METHODS 

In this section we present the general procedure for fitting model 
light-curves with empirical templates and for inferring intrinsic 
photometric properties from observed SN la light curves. We also 
present the algorithm used to cross-correlate the synthetic spectra 
with a large database of observed SN spectra, which we use as a 
basis for quantitatively evaluating each model. 



3.1 Light-curve fits 

We deri ve synthetic UBVRI (■Bessell"l990^ and JHKs (2MASS 
system; ICohen, Wheaton & Megeath 2003) magnitudes from the 
model spectra, with associated errors based on the number of pho- 
ton packets in a given frequency bin. Typical errors at maximum 
light are < 0.005 mag in UBVR, < 0.01 mag in /, and pro- 
gressively larger errors (0.015-0.05 mag) for JHKs. We show the 
UBVRIJHKs fight curves for model DD2D_iso_06_dc2 viewed 
along ^ = 88 ° in Fig. [ 3 conipared to observations of SN 2003du 
dAnupama. Safiu & Jos?2005l : lLeonard et alJ2005l:IStanishey et alJ 



l2007h and the NIR template fight curves of iMandel et"aD (l2009h 
The models appear to match the data well in the optical (albeit with 
significant differences in the U and / bands), but deviate in the near 
infrared. 



3.1.1 Direct polynomial fits 

We determine the time of maximum light and peak magnitude in 
the optical bands (UBVR), as well as the corresponding Amis de- 
cline rate. Our final estimates for the various quantities, along with 
their associated errors, are based on 1000 Monte Carlo realizations 
using the magnitude errors associated with each light-curve point. 
The resulting error on the fit parameters is below 0.1% in all cases. 



3.1.2 Template-based fits 

We attempt to fit BVRI light curves for each model viewed from 
all 30 directions with three different light-cu rve "fitters" commonly 
used amongst SN observers: MLCS2k2 (Jha, Riess & Kirshner 
l2007h . SALT2 dOuv et al.ll2007h . and SNooPy (Burns et al. 2011). 
To first order, these algorithms all share the common approach of 
comparing an input light curve with empirical templates based on 
observed data. Adding the [/-band light curves for the MLCS2k2 
and SALT2 fits (there are no standard [/-band templates for 
SNooPy) has negligible impact on the results. Since no extinction is 
applied to the model synthetic magnitudes, we force the extinction 
to zero. This is not possible with S ALT2, as the effects of extinction 
and intrinsic colour variations are described by a single parameter. 

We show an example MLCS2k2 fit to the B-band light curve 
for model DD2D_iso_06_dc2 viewed along = 46° in the left 
panel of Fig. |4] (blue hatched area). The fit is formally accept- 
able overall, with a reduced xl — 0.79. However, a closer inspec- 
tion shows that the MLCS2k2 template reaches its peak magnitude 
'^ 2.5 d earlier than the model predictions. This suggests that the 
models rise more slowly to maximum light than actually observed. 
In MLCS2k2 we have the possibility of fixing the time of maxi- 
mum in the B band, and we use the value found from direct poly- 
nomial fits to the model light curves. The fits are severely degraded 
(xl — 3.37; orange hatched area), and the empirical templates 
again appear to require a faster rise time than the models suggest. 

This tension between model and empirical light-curve shapes 



is not specific to MLCS2k2, as can be seen from the compari- 
son with SALT2 and SNooPy on the B-band light curve for the 
same model/viewing angle combination (Fig. IH right panel). In 
this particular example, the time of B-band maximum is underes- 
timated by approximately 3 d, 2 d, and 1 d for MLCS2k2, SALT2, 
and SNooPY, respectively. The fits become formally worse follow- 
ing this same sequence, but this is due to the progressively smaller 
errors associated with the empirical templates (the mean absolute 
deviation of the template from the synthetic light curve over the 
time interval -15 ^ t [d] ^ +20 is ^ 0.7 mag for the three fit- 
ters). 

A comparison of light-curve fits for all model/viewing angle 
combinations shows that, on average, MLCS2k2 underestimates 
tmax(5) by - 2d, SALT2 by - 0.5 d, and SNooPy by - 0.3d. 
One might then be tempted to assume that the models have an over- 
all shape that is compatible with the SNooPy empirical templates, 
but a comparison of Amis (B) values inferred from these fits with 
those directly measured on the model light curves reveals a strong 
systematic trend (despite the residuals having a mean consistent 
with zero; see Fig.[5]). The result is the same when we fix tmax{B) 
to its actual value. 

Given the difficulty for empirical templates to correctly re- 
produce the overall shape of the model light curves, in what fol- 
lows we use the quantities derived from the direct polynomial fits 
(time of maximum, peak magnitude, post-maximum decline rate) 
for the models. The observed light curves are not always suffi- 
ciently well sampled to enable us to determine the time of max- 
imum and Amis(i?) directly from polynomial fits, in which case 
we use values inferred from light-curve fitters. For SNe la for which 
we have both direct polynomial and template-based fits, the differ- 
ence in time of maximum and decline rate from both methods is 
consistent with zero. 

The various observables derived from the models are in- 
trinsic to the supernova. Observations, however, are affected by 
extinction, both in the Milky Way and the SN host galaxy. 
The former is usually determined using the dust maps of 
ISchlegeL Finkbeiner & DavisI (Il998l) , but the latter requires as- 
sumptions about the dust properties and the intrinsic colours of 
SNe la. In the next section we describe the BayeSN statistica l 
model for SN la light curves (iMandeL Naravan & Kirshnej|20Ilh 
used to infer the host-galaxy extinction and hence the intrinsic 
SN la peak magnitudes and colours. 



3.1.3 Inferring intrinsic SN la properties with BayeSN 

iMandeL Naravan & KirshneJ (l201lh constructed hierarchical 
Bayesian models for SN la light curves spanning optical through 
near infrared (BVRIJH) data. These statistical models describe 
the apparent distribution of light curves as a convolution of 
intrinsic SN la variations and a dust distribution. They modeled 
the intrinsic covariance structure of the full multi-band light 
curves, capturing population correlations between the intrinsic 
absolute magnitudes, intrinsic colours and light curve decline 
rates over multiple phases in the optical and NIR filters, as well 
as the distribution of host galaxy dust and an apparent correlation 
between the dust extinction Ay and its wavelength dependence, 
parameterized by Ry. The models use the optical and NIR time 
series data to fit individual SN la light curves, estimate their dust 
extinction and predict their distances. Probabilistic inference of 
the parameters of individual SNe la and also those describing the 
intrinsic SN la and dust populations are computed using a novel 
Markov Chain Monte Carlo code. These models were trained on a 
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nearb y {z < 0.07) set of SNe l a with optical (CfA3,rHicken et all 
I2OO9I : Carn epe SN Program JContreras et aP l201Qh and NIR 
fPAIRITEL.' Wood-Vasev et al.l l2008') data, plus light curves from 
the literature with joint optical and NIR observations. The resulting 
Markov chain is used to estimate the dust extinction, apparent and 
absolute light curves, and intrinsic colours for each SN la. For this 
study, we employed the model that captures a linear correlation 
between Ry^ and A y, as this produced the best c ross- validated 
distance predictions in lMandel. Naravan & Kirshnej (12011.) . 



of the template spectra (see ^ 12.21 ). Last, we restrict the ages of the 
template spectra to be within 3 d from B-band maximum light (re- 
laxing this age constraint does not significantly alter the results at 
maximum light). We are then left with 441 spectra of 175 SNe la 
in the SNID database. 

The strength of a correlation between the input model spec- 
trum and a particular template is embodied in the rlap parame- 
ter, which is the product of th e height of the norm alized cross- 
correlation peak (the r- value of Tonry & Davis 1979) and the over- 
lap in log- wavelength space between the input and template spec- 
tra (here lap ^ ln(7500/3500 ) ^ 0.7). A correlation is consid- 
ered good when rlap ^ 5 fsee lBlondin & Tonrvll2Q07n . For each 
model/viewing angle combination, we record the best-match tem- 
plate (the one with the highest rlap), as well as the mean rlap value 
for the top five matches, which we consider a more reliable esti- 
mate of how similar a particular model spectrum is to observations. 
Fig.[6lshows some example SNID fits to our subset of eight models 
(see Table [B, viewed along the direction which provides the best 
match to an observed spectrum. The matches are good overall, with 
some discrepancies around 4000 A for some models. 

We have also run SNID on model spectra at ages other than 
maximum light, namely — 10 ^ t [d] ^ +20, each time impos- 
ing a ±3 d age constraint on the template spectra. Close to maxi- 
mum light (— 5 ^ t [d] ^ +5), the fits are of similar quality as at 
maximum, i.e. the mean rlap value for the cross-correlation with 
observed spectra is on average the same. However, further away 
from maximum light, the fits are degraded, and in some cases no 
good matches are found. Relaxing the age constraint for the tem- 
plate spectra results in formally better matches (higher mean rlap), 
but the mean age of the top five matches tends to over- (under-) esti- 
mate the actual age for model spectra before (after) maximum light. 
This again suggests some differences between the rate at which the 
synthetic spectra evolve compared to observed spectra. 



3.2 Spectral comparison using SNID 

We use the S uperN ova IDentification (SNID) code of 
iBlondin & Tonrvl (120071) to cross-correlate the synthetic spec- 
tra with a large database of observed spectra (referred to as 
"templates"). SNID is commonly used to determine the type, 
redshift, and age of a supernova spectrum (e.g. in lAU circulars). 
The database associated with the public release of SNID includes 
493 spectra of 48 SNe la in the age range -15 ^ t [d] ^ +30. 
We have augmented it with new data from the CfA SN Program 
(see ^ [^21 and r e cently published SN la data (SN 2 005bl : 
Taubenberger et alJ I2OO8I : SN 2005cf: iGaravini et alJ l2007l . 



Wang et al. 2009b ; SN 2005cg: lOuimbv et al.ll2006l:"sN 2005hi: 
Ouimbv. Hoflich & Wheeled l200i SN 2007ax: iKasliwal et all 



2008h . 



This revised SNID database now includes 2046 spectra 
of 274 SNe la o ver the same age range. We refer the reader to 
iBlondin & Tonrvl (12007 ) for an extensive discussion of the SNID 
algorithm. We point out here that both the input and template 
spectra are "flattened" through division by a pseudo continuum, 
such that the correlation relies on the relative shape and strength of 
spectral features, and is therefore insensitive to colour information, 
which can be intrinsic, or result from extinction by dust and 
flux-calibration uncertainties. 

We first run SNID on the model spectra closest to B-band 
maximum light. We only consider SN la spectral templates, and 
force the redshift (usually a free parameter in SNID) io z — 0. We 
restrict the rest-frame wavelength interval over which the correla- 
tion is done to 3500-7500 A to match the wavelength range of most 



3.3 Model selection and the width-luminosity relation 

Based on the SNID fits to the maximum-light model spectra, we 
can assign a "grade", namely the mean rlap value for the top five 
matches, to each model/viewing angle combination. Fig. |7] {left) 
shows the width-luminosity relation (Mmax(i^) vs. Ami5{B)) for 
the 2D delayed-detonation models of KRW09 (see their Fig. 3), 
where each point has now been colour-coded according to the 
mean rlap value for that particular model/viewing angle. Points 
with higher rlap values tend to lie systematically on the observed 
width-luminosity relation, i.e. the agreement between synthetic and 
observed spectra appears to translate into an agreement in photo- 
metric properties, thereby justifying our ranking scheme based on 
maximum- light spectra only (but see ^ 14.31 ). 

Since each model has equal probability of being observed 
from one of the 30 different viewing angles (see § [3, it can only 
be validated if all viewing angles yield a good match to an ob- 
served SN la spectrum. Conversely, a lack of a good match for 
a single viewing angle should in principle lead us to reject the 
model as a whole, by which we mean it cannot be considered a 
valid approximation of observed SNe la. Six models fall in this cat- 
egory: DD2D_asym_03 (dc2 and dc3), DD2D_asym_06 (dc2 and 
dc3), and DD2D_asym_08 (dc2 and dc3). As their names sug- 
gest, all are models with an anisotropic distribution of ignition 
points. Moreover, they all have slow-declining light curves, with 
0.76 ^ Ami5{B) ^ 0.89. While the DD2D_asym_03 models (dc2 
and dc3) have 1-2 viewing angles with no good SNID matches, both 
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DD2D_asym_06 and DD2D_asym_08 each have no good matches 
for more than half of the viewing angles. 

We have selected a subset of eight models for detailed stud- 
ies (Fig. [T] right): DD2D_iso_03_dcl (red) is the model with the 
highest rlap value overall; DD2D_asym_01_dc3 (orange) is the 
best model with a highly asymmetric distribution of ^^Ni (^ 70% 
variation in peak bolometric luminosity; see Table [B; the com- 
plete DD2D_iso_06 model series (del through dc5) has consis- 
tently high rlap values and spans a large range of Ami5{B) 
(0.8-1.6 mag), thus enabling us to isolate the impact of the crite- 
rion for deflagration-to-detonation transition. It also includes the 
reference model DD2D_iso_06_dc2 used by KRW09 for detailed 
post-processing nucleosynthesis calculations. These seven models 



all lie on the observed WLR and span a large range in synthe- 
sized ^^Ni mass (0.42-0.88 Mq). Last, DD2D_iso_08_dc3 (green; 
lower right) is the best model that does not lie on the width- 
luminosity relation. For all models other than DD2D_iso_03_dcl 
and DD2D_iso_08_dc3, ^ 15° corresponds to the direction of 
largest extent of ^^Ni distribution, and hence to the most luminous 
point in each model series. This is in part a consequence of the sym- 
metry imposed on the explosion by the 2D axisymmetric setup in 
the models of KRW09. These models will be highlighted and their 
colour-coding preserved in subsequent figures. 

In what follows we compare photometric (§ |4]) and spectro- 
scopic (§|5]) properties separately, as is commonly done in observa- 
tional studies of SNe la. This in part reflects the way the observa- 
tions are done (generally there are separate follow-up programmes 
for photometric and spectroscopic observations) and calibrated (in- 
dividual photometric measurements have a typical absolute cali- 
bration accurate to a few percent, whereas individual spectra have 
a relative flux calibration that is generally less accurate; see, e.g., 
[Matheson et al.l 12008). In principle one would want to compare 
light curves and spectra simultaneously, and in fact we will refer 
to conclusions drawn in §|4]when discussing spectroscopic proper- 
ties in § [5] 



4 COMPARISON OF PHOTOMETRIC PROPERTIES 

In this section we focus on multi-band light curves only, by com- 
paring the rise times, maximum-light colours and colour evolution 
in models and data. 



As noted earlier (§ 13.11 ), the models rise too slowly to maximum 
light in the B band. We show the relation between the B-band 
rise time (trise{B)) and AmiB{B) in Fig.[8l for the 2D models of 
KRW09 (upper panel) and as measured on SN la light-curves from 
the SDSS-II Supernova Survey bv lHavden et all (12010) (where we 
exclude SNe la with Ami5(B) > 1.6 to match the Ami5(B) 
range of the models). In both cases trise{B) appears to be largely 
independent of Ami^{B). The range of model rise times (18- 
23 d) overlaps with the observed distribution (13-23 d), although 
the bulk of the data display rise times smaller than 18 d. Within a 
given model, the rise time changes with viewing angle according 
to the distribution of ^^Ni along that particular inclination, any- 
where between ^ 1 d and ^ 3 d. For model DD2D_iso_06_dc2, the 
B-band rise time is ^ 20 d regardless of viewing angle, while for 
model DD2D_asym_01_dc3, it varies between - 20d (6* < 90°) 
and -23 d (6* = 165°). 

The model B-band rise times are more consistent with earlier 

II 1 

Riess et al. 1999a), but 
the latter were based on empirical templates that were "stretched" 
by equal amounts in the rising (pre-maximum) and fa lling (post- 
maximum) portions of the light curve. Istrovinki (120071) considered 
asymmetric templates and found a mean rise tim e of 17.4 zb 0.4 d, 
subsequently confirmed bv lHavden et al.l (I2OIOI) using a larger data 
set. 

The variation in rise times between various bands is linked 
to the redistribution of flux in frequency in the supernova ejecta. 
Fig. |9] shows the difference in time of maximum light with respect 
to the B band, for the UVRIJ photometric bands as well as for the 
integrated UBVRI (UVOIR) flux, where we have used polyno- 
mial fits to observed SN la UBVRIJ light curves (again exclud- 
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Figure 8. B-band rise time v^. Ami5(B), for the 2D delayed-detonation 
models of KRW09 (top) and for SN la data from Havden et al. (2010) (bot- 
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subset of selected models (filled circles). The open circles correspond to 
all other models. The rotated histogram in the lower panel shows the dis- 
tribution of the Hay den et al. measurements, while the vertical error bar 
corresponds to the range of model rise times. 



ing SNe la with Ami5{B) > 1.6) from the literature for the data 
distributions (open histograms). The /-band and NIR light curves 
of SNe la display characteristic secondary maxima; here we report 
the time corresponding to ihQ first maximum. The times of UVOIR 
maximum were derived from light curves published by .Stritzinger. 
OOOSh . where we used only SNe la with UBVRI light curves so 



as to minimize the uncertainty introduced by bolometric correc- 
tions which account for missing [/-band data. These distributions 
ar e broadly consistent with those found o n a more limited data set 
bv lContardo. Leibundgut & Vaccal feOOOh (their Fig. 4). 

The data indicate that SNe la peak i n the U b and typically 
^ 2 d earlier than in B (consistent with Uha et al.l (12006.) . who 
found tnaax([/) - tnaax(5) = -2.3 =b 0.4 d), and - 1-2 d later 
in V and R. The /-band has a large spread in times of (first) max- 
imum, but most SNe la peak ^ 2d earlier than in B. The redder 
NIR bands peak even earlier (only J is shown in Fig. [9]), as pre - 
viously noted by several authors (e.g. jMeikle & Hernandezll200Qh . 
Where a simplistic model would predict the reddest bands to peak 
the latest as the temperature of the ejecta decreases (and this ex- 
plains the tmax([/) < tmax(i^) < tmax(V) scqucnce), thc fact that 
the IJHKs bands peak before the B band shows that flux is effi- 
ciently redistributed from the ultraviolet to the near infrared (e.g., 
IPinto & Eastmanll2000bh . The magnitude of this effect is strongly 
dependent on the ionization stage of iron-group elements in the 
SN la core, and is responsible for the secondary maxima observed 
in IJHKs SN la light curves rsee lKasenll200d) . 

This trend is broadly reproduced by the models, with a large 
overlap with the observations in VRIJ, and a more marginal over- 
lap in U (especially when considering our subset of selected mod- 
els). The distribution of the difference between UVOIR and B-band 
maxima in the models, however, has no overlap with the observed 
distribution. The data seem to indic ate that the times of UVOIR and 
B-band maxima coincide (see also lContardo. Leibundgut & Vaccal 
2000), whereas the model bolometric light curves peak 1-2 d ear- 
lier than the B-band. The U band contributes most to this (modest) 
discrepancy, as it is the brightest band at maximum light in these 
models. 



4.2 Maximum-light colours 

Fig. [TOl shows the normalized distributions of intrinsic U — B, 
B -V, B - R,V - R,V - I^andV - J colours at B-band maxi- 
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mum as inferred from BayeSN fits to actual data {open histograms) 
compared with the ranges for the 2D models of KRW09 {gray, all 
models; black: subset of selected models). The data indicate that 
SNe la have an intrinsic B — V ^ —0.1 mag at B-band maximum, 
whereas the models have a systematically redder B — V colour 
(0.0-0.2 mag). This was already noted by KRW09 and attributed 
to the approximate treatment of non-LTE effects. The distributions 
fox B — R,V — R, and V — I have a large overlap between the 
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bars). The U — B colour distribution is derived from MLCS2k2 fits to 
SNe la with low visual extinction {Ay < 0.1 mag). 

models and data, but contrary to ih^ B — V colour the former have 
a much larger spread (this is also true when only considering the 
subset of selected models). The V — J colour displays differences 
of up to '^ 0.5 mag between the models and the data, albeit with a 
large overlap with the observed distribution. 

The models have a systematically redder B — V colour at B- 
band maxi mum, but they are systematically blu er inU — B colour 
(Fig. [To]). MandeL Narava n & Kirshnej (1201 ih did not include U- 
band data in their analysis to avoid calibration uncertainties, so we 
instead infer the intrinsic [/-band magnitude based on MLCS2k2 
fits to SNe la with low visual extinction {Ay < 0.1 mag). The mod- 
els show some overlap inU — B colour with the data, but extend 
> 0.2 mag further in the blue. This excess of [/-band flux can be 
explained either by a hotter ejecta temperature affecting the overall 
SED or an incomplete description of line-blanketing of UV pho- 
tons and subsequent redistribution of flux to redder wavelengths, 
or a combination of both. Since the impact of such effects is most 
easily seen on the spectra, we will discuss them in §|5] 

More luminous SNe la are also int rinsically bl uer (this is 
known as the "brighter-bluer" relation, e.g. lTripplll998h . Its degen- 
eracy with extinction by dust (SNe la obscured by dust will appear 
redder) is exploited by some light-curve fitters to describe the ef- 
fects of extinction and intrinsic colour variations using a single pa- 
rameter (e.g., SALT2: Guv et al. 2007) . The brighter-bluer relation 
(i.e. Mmax(i^) vs. intrinsic B — V ai B-hmid maximum) derived 
from BayeSN fits to actual data is shown in the upper-left panel 
of Fig.\ll\{dashed line). The 2D models of KRW09 also follow a 
brighter-bluer relation, but it is significantly steeper than observed. 
This difference is not surprising given the lack of overlap in B — V 
colour between the models and the data (Fig.fTol). 

The agreement between models and data is significantly bet- 
ter when considering Mma^{B) vs. B — R at B-band maximum 
(Fig. [TT] upper-right panel), which forms an extremely tight rela- 
tion. Small deviations are apparent forB — R > Omag, where 
the models predict peak colours that are much redder than inferred 
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data used to determine the observed relation, while the grey hatched area corresponds to extrapolations of this relation. 



from observations (Fig. [TO]). Of all the selected models, model 
DD2D_asym_01_dc3 (filled orange circles) displays the largest 
variation of B — R colour (the same is true to a lesser extent when 
considering the B — V colour), illustrating the greater change of 
the overall SED with viewing angle stemming from the asymmet- 
ric ^^Ni distribution (see Fig.[T]). 

The lower panels of Fig. [TT] show Mr and Mi vs. intrinsic 
R — I colour (all quantities are at B-band maximum). While Mr 
has almost no relation io R — I colour (the slope is consistent with 
zero at ^ la), Mj vs. R — I appears to define a bnghter-redder 
relatioru. In both cases, models that predict redder R — I colours 
than inferred from observations lie off the relation. 



4.3 Colour evolution 

The early-time evolution of SN la colours results from a complex 
interplay between a gradual cooling of the ejecta and changes in 
ionization stages of several species, leading to order-of-magnitude 
changes in the opacity at selected frequencies. In particular, the 
III^II recombination of iron-group elements is responsible for 
most of the redistribution of flux from blue to red wavelengths and 
its timesc ale controls the rate of colour evolution in SNe la (see 
iKasen & Wooslev__2007_\ 

^ While this is not the subject of this paper, such a relation would help 
break the degeneracy between extinction by dust and intrinsic colour varia- 
tions. 



Fig. [12] shows the time evolution of U - B, B - V, B - R, 

V — R,V — I and V — J colours for our subset of selected models 
(we have chosen one representative model for the DD2D_iso_06 
series, namely DD2D_iso_06_dc2), colour-coded according to the 
viewing angle 0. Overplotted are observed (dereddened) colour 
curves from SNe la in similar AmiBJB ) ranges, except for the 

V — J colour curves where we show the iKrisciunas et al.! (120041) 
templates for slow- (0.8 < /\mi^{B) < 1.0) and mid-range 
(1.0 < AmisiB) < 1.4) decliners, as well as for SN 2001el 
(Ami5{B) ^ 1.1). To minimize the uncertainties associated with 
extinction corrections, we only include SNe la with a visual extinc- 
tion Av < 0.25 mag, as inferred from BayeSN fits. 

The B — V colour evolution of the models is largely con- 
sistent with that seen in the data, with small systematic offsets 
to redder colours for model DD2D_iso_08_dc3, which lies off the 
width-luminosity relation. The agreement between the models and 
the data for the B — R colour evolution is also satisfactory, apart 
for model DD2D_iso_03_dcl, whose B — R colour becomes pro- 
gressively bluer with age until ^ 10 d past B-band maximum while 
the data become redder over the same time interval. The disagree- 
ment is exacerbated when considering V — R for this model, but 
the V — I colour evolution is in good agreement (as is true for 
the other models of Fig. [12]). Interestingly, DD2D_iso_03_dcl is 
our best overall model based on cross-correlations with observed 
spectra using SNID (see § 13. St . Since SNID compares the rela- 
tive shapes and strengths of spectral features, the overall SED was 
not taken into account when ranking these models, which allows 
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Figure 12. Time evolution of intrinsic colours measured at B-band maximum for four of our selected models, colour-coded according to the viewing angle 
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for possibly large offsets in broadband colours between our se- 
lected models and the data. Based on iheV — R colour evolution, 
one would rank model DD2D_iso_08_dc3 highest of the four mod- 
els shown in Fig. [TJl where in fact it is the only model that does 
not lie on the observed width-luminosity relation. As noted earlier, 
model DD2D_asym_01_dc3 displays the largest colour variation at 
any given time. 

The colour curves forU — B and V — J, however, reveal the 
largest inconsistencies between the models and the data. All the 
models have a systematically bluer U — B colour at early times 
(up until ^ 10 d past B-band maximum), and increasingly redder 
U — B colours at later times. By 20 d past B-band maximum, the 
models display Si U — B colour that is ^ 0.5 mag too red with 



respect to the data, a difference that clearly cannot be explained by 
simple [/-band calibration uncertainties. 

The qualitative behaviour of the F — J colour curves is 
compatible with the data, but displays large offsets (> Imag 
for model DD2D_iso_03_dcl) around ^ lOd past B-band maxi- 
mum, which corresponds to the time of flux minimum between the 
two J-band maxima (see Fig. [3]). This is not the case for model 
DD2D_iso_08_dc3, which has redder V — J colours at any given 
time when compared to observations. The failure of the models to 
properly reproduce the NIR flux in SNe la at all times could be 
due in part to the atomic data used in the radiative transfer cal- 
culations, as noted bv iKasenI (l2006h . By comparing J-band light 
curves computed with the 500,000 line Kurucz CD 23 and the 42 
million line Kurucz CD 1, he showed that including more lines re- 
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suited in '-^ 1 mag more J-band flux between —10 and +20 days 
from B-band maximum and a shallower flux minimum (his Fig. 5), 
although it is not clear which lines contribute most to this effect. 
Other approximations used in the radiative transfer calculations 
(LTE, expansion opacity formalism, equivalent two-level atom for 
flux redistribution etc.) could also explain part of these discrepan- 
cies. 

Fig. [12] conveys the message that a direct comparison 
of SN la models with observations must involve both photo- 
metric and spectroscopic properties. Our highest-ranked model 
based on cross-correlations with spectral templates using SNID 
is DD2D_iso_03_dcl, yet its colour curves present the largest 
level of disagreement with observations. We also note that 
DD2D_iso_08_dc3 fares best based on the evolution of its broad- 
band colours, yet it lies off the observed width-luminosity relation. 



5 COMPARISON OF SPECTROSCOPIC PROPERTIES 

It is difficult to resolve discrepancies between the models and data 
based on photometric measurements alone, and we now turn to the 
comparison of spectroscopic properties, first considering spectra at 
maximum light (^ 15. IK and then their evolution with time (^ 15.21 ). 



5.1 Maximum-light spectra 

5.1.1 Spectroscopic diversity at a given /\mi^{B) 

We would like to compare the general appearance of the model 
spectra in specific Amis (5) ranges. To do this we select two 
narrow Ami^iB) ranges (±0.05 around a central value) that in- 
clude at least one viewing angle for each of the eight models in 
our subset. Based on Fig.[71 we see that three of our subset of se- 
lected models overlap at Ami5{B) ^ 0.95 (DD2D_iso_03_dcl, 
DD2D_iso_06_dcl, and DD2D_asym_01_dc3), while the remain- 
ing five overlap at Am 15(B) ^ 1.3 (DD2D_iso_06_dc2 through 
dc5 and DD2D_iso_08_dc3). We thus define our Ami5{B) ranges 
of interest as 0.90-1.00 and 1.25-1.35. Figure [13] shows observed 
maximum-light spectra in both Ami^iB) ranges (black) compared 
to our subset of selected models (red), for which we only show 
spectra for one viewing angle (the one that yields a Amis (B) clos- 
est to the middle of the Amis(B) range considered). The observed 
spectra all have an estimated visual extinction Ay < 0.25 mag and 
have been dereddened. 

For 0.90 ^ Amis(B) ^ 1.00, the observed sample con- 
sists of both "normal" and 1991T- like SNe la (of which SN 1998es, 
SN 1999aa, SN 1999dq; all from lMatheson et al.ll2008h . The latter 
are characterized by hot ejecta, with evidence for lines of doubly- 
ionized iron-group elements (Fe Ill/Co III) in place of lines of 
Fe Il/Co II in maximum-light spectra of "normal" SNe la, and 
which explains part of the observed variation around ^ 4300 A. 
Redward of this, the spectra are remarkably similar to one another, 
albeit with variations in the strength of the absorption component 
of the prominent Si II 16355 line (1991T-like SNe la have shallower 
Si II absorptions). The model spectra in the same Amis(i?) range 
are also similar to one another, but they have a much bluer (hotter) 
SED than the data (as noted based on the comparison ofU — B 
colours in ^ I4.2K and shallower Si II 16355 absorptions. This line 
is also broader and its absorption component more blueshifted than 
observed in the data, both evidence for larger ejecta velocities. Sev- 
eral small-scale features present in the data are thus "washed out" 
in the models, as the higher expansion velocities increase the line 



overlap. This is clearly seen in the S II 115454,5640 doublet, which 
appears as two distinct absorption components in all the observa- 
tions but which is barely resolved in the models (and not resolved 
at all for model DD2D_iso_03_dcl). 

For 1.25 ^ Amis(B) ^ 1.35, the observed sample con- 
sists almost exclusively of "normal" SNe la, apart from SN 2004dt 
which displays larger absorption v elocities (cl e arly v isible in the 
Si II 16355 line) and was found bv lWangetaP (l2006h to have the 
highest degree of polarization ever measured in aSN Ia. SN 2004eo 
was labeled "transitional" by Pastorello e t al] (12007!) due to its in- 
termediate properties between normal and underluminous SNe la. 
The model spectra in the same Amis(B) range show a large de- 
gree of heterogeneity, both in the slope of the overall SED and 
in the relative shapes and strengths of spectral features. Model 
DD2D_iso_06_dc2 is too blue compared to observations, and dis- 
plays relatively shallow absorption features; conversely, models 
DD2D_iso_06_dc4 and dc5 have a significantly redder SED with 
broader and deeper absorption features than in the observations, as 
clearly seen from the three prominent Si II absorption features at 
- 4000 A, - 5800 A, and - 6100 A. Model DD2D_iso_08_dc3 
fares better both in terms of SED slope and relative strengths 
of spectral features, albeit with less prominent absorption around 
^ 4300 A and lack of small-scale structure around ^ 4800 A (both 
features are predominantly due to lines of Fe Il/Co II, with im- 
portant contributions from Mg II and Ti II). As was the case for 
the model spectra with smaller Amis (5) (Fig. [13] left panel), the 
S II 115454,5640 doublet is not resolved in this model (the same is 
also true for DD2D_iso_06_dc2). 

The partial failure of the models to reproduce the relative 
shapes and strength of the Si II lines impacts their ability to re- 
produce the observed correlati on of several spectroscopic indica - 
tors with Amis(B) (see, e.g., Blondin. Mandel & K irshner 2011). 
Fig. [T4l shows two such indicators, namely the 7^ (Si) ratio of 
iNugent et al.l (Il995n (defined as the ratio of the relative absorption 
depth of the Si II 15972 line to that of Si II 1635 5) and the pseudo 
equiv al ent width (pEW) o f Si II 14130 (see lArseniieyic et alj 
2008; Walker etal. 2010; B londin. Ma ndel & Kir shnej 1201 ll : 
Chotard et al. 2011). There is a clear correlation between both in- 
dicators and Ami5(B) in our spectroscopic data set (lower pan- 
els), albeit with two significant outliers for 7^ (Si) (one of which 
SN 2006bt, whos e peculiar nature has been discussed at length by 
Folev et al]|20ld ). A linear fit to the data is shown in both cases 
(dashed line) and overplotted on measurements from the models 
(upper panels). There is no clear correlation between 7^(Si) and 
Amis (B) in the models, but the measurements have a large associ- 
ated error due to Monte Carlo noise affecting the weaker Si II 15972 
line. As for Si II 14130, there is a general trend of increasing pEW 
with Amis(B), but with a large scatter for Amis (5) ^ 1.1 
caused by some model spectra with unusually broad and deep ab- 
sorptions. 



5.1.2 Composite spectra at a given Attiis (B) 

To better study overall systematic differences between the models 
and observations, we generate composite spectra using the same 
pre-processing as done by SNID: the individual spectra are "flat- 
tened" through division by a pseudo continuum. We then compute 
the mean flux in each wavelength bin, as well as the standard devi- 
ation from the mean. The result is a composite spectrum with error 
bands, which we show in Fig. [15] Note that the composite spectra 
based on observations (hatched grey) only uses the sample of spec- 
tra shown in Fig. [13] while those based on the models comprise all 
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Figure 13. Maximum-light spectra for 0.90 ^ /S.mi^{B) ^ 1.00 (left) and 1.25 ^ /S.mi^{B) ^ 1.35 (right) compared to models in the same /S.mi^{B) 
range. For sake of clarity, we only show spectra corresponding to our subset of selected models, and select the viewing angle corresponding to the decline 
rate closest to the middle of the Ami5(B) range considered. The vertical lines correspond to S ii 15640 {dotted) and Si ii 16355 {dashed) blueshifted by 
10000 km s-^ 



viewing angles which yield a Ami5{B) in the appropriate range. 
The composite spectra for all (selected) models thus consists of 234 
(50) and 150 (41) individual spectra for 0.90 ^ Amis (5) ^ 1.00 
and 1.25 ^ Ami5{B) ^ 1.35, respectively. 

The composite spectra for 0.90 ^ Amis (5) ^ 1-00 confirm 
the shallower absorptions and larger blueshifts of spectral lines in 
the models. The location of maximum absorption in the Si II 16355 
line is typically ^ 100 A too blue in the models, corresponding to a 
^ 5000 km s~^ difference in absorption velocity. The small-scale 
structure in the two iron-group-dominated absorption features (^ 
4300 A and ^ 4800 A) present in the data is not seen in the models, 
as mentioned earlier. The emission component of the Si II 16355 
line appears to be relatively stronger in the models, but this is an 
artefact of the division by a pseudo continuum. The flux level in the 
models drops significantly redward of ^ 6500 A compared to the 
data, artificially enhancing the relative strength of the Si II 16355 
emission profile after division by this pseudo continuum. 

The composite spectra for 1.25 ^ Amis(B) ^ 1.35 do not 
seem to support the apparent heterogeneity of the model spectra in 
Fig. [13] This is not surprising since Fig. [13] only shows 5 of the 



150 (41) spectra used to generate the composite spectrum for all 
(selected) models. It is also possible that the division by a pseudo 
continuum attenuates some of these differences. Nonetheless, the 
composite spectra for this Amis (5) range do reveal some interest- 
ing properties of the models. The absorption velocity offset is less 
pronounced in the Si II lines. The location of maximum absorption 
in the Si II 14130 line even appears to be at larger wavelengths in the 
models, but this could result from line overlap. The two iron-group- 
dominated features at ^ 4300 A and ~ 4800 A again appear to be 
smoothed by a larger expansion velocity field. The S II doublet is 
again shallower in the models, although better resolved than for 
0.90 ^ Am 15(B) ^ 1.00. Both absorption components appear to 
be slightly offset to redder wavelengths when compared to the data, 
although the difference is typically < 25 A, or < 1500 km s~^. 
The emission component of the Si II 16355 line is again artificially 
enhanced via the division by a pseudo continuum, although part of 
the difference is real and clearly seen in the model spectra shown 
in Fig.[T3](in particular for models DD2D_iso_06_dc4 and dc5). We 
will explore some of these systematic differences in more detail by 
focusing on this line in the next section. 
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Figure 14. Left: Spectroscopic indicator 7^(Si) v^. Ami5(B), for the 2D delayed-detonation models of KRW09 (top) and for SN la data from the CfA SN 
Program (bottom). The colour-coding is the same as in the right panel of Fig. [7] The dashed line is a linear fit to the CfA data, and is overplotted in the upper 
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Figure 15. Composite maximum-light spectra for 0.90 ^ /S.mi^{B) ^ 1.00 (left) and 1.25 ^ /S.mi^{B) ^ 1.35 (right). The shaded bands correspond to 
the standard deviation about the mean maximum-light spectrum. We show the flattened composite spectra for all DD2D models (filled dark blue), our subset 
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5.1.3 Si II t6355 absorption velocity 



IWang et alJ (12009a) introduced a classification scheme for SNe la 
based on the Si II 16355 absorption velocity at maximum light. 
By considering deviations from an empirical mean trend, they 
classified their SN la sample into "normal" and "high-velocity" 
subclasses, with separate subclasses for 1991T- and 1991bg-like 
SNe la. Fig. [161 (lower panel) shows the relation between the 
pseudo equivalent width of the Si II 16355 line with its absorption 
velocity (i.e. the velocity at maximum absorption, ^;abs) within 3 d 
from maximum light in our spectroscopic sample from the CfA SN 
Program (see Wang et al. 2009a, their Fig. 2). Only SNe la with 
Am 15 (5) < 1.6 are included, to match the Am 15(B) range of the 
models. We added a random component drawn from a 500 km s~^- 
wide uniform distribution to account for velocity measurement er- 



rors resulting from the 10 A binning in the synthetic spectra. We 
have also run one model (DD2D_iso_06_dc2) with a larger num- 
ber of photon packets (10^ instead of 10^) and find the measure- 
ments to be consistent with the lower-S/N model run. Almost all 
the models have Si II 16355 absorption velocities consistent with 
the "high-velocity" subclass, regardless of the viewing angle, in 
agreement with the large blueshifts discussed in the previous sec- 
tion. Several models have significantly lower pEW than in the data 
(at a given i^abs), again consistent with the shallower absorptions 
seen in some of the individual spectra (in particular, those in the 
0.90 ^ Ami5(B) ^ 1.00 range). 

These measurements reveal further properties of the mod- 
els themselves. Model DD2D_asym_01_dc3 (filled orange circles) 
displays the largest variation in Si II 16355 absorption velocity 
(-11000 > i;abs > -15000 km s"^), the larger blueshifts cor- 
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responding to the largest extent of the ^^Ni distribution. For mod- 
els with a symmetric distribution of ^^Ni, one expects the expan- 
sion rate to scale approximately with ^/E^[n/M. Fig. [16] shows 
that one of our selected models with the lowest kinetic energy 
(DD2DAsoJ)6.dc5, filled yellow circles; ^kin ?^ 1.3 x 10^^ erg, 
see Table [T]) has lower |i;abs| than one of our selected models with 
the highest kinetic energy (DD2D_iso_03_dcl, filled red circles; 
Emu ^ 1.5x 10^^ erg). To see whether this holds for the other 2D 
models of KRW09, we show the relation of ^;abs,max(Si II 16355) 
vs. kinetic energy in Fig.fTvl 

We do not see a general trend of increasing |i;abs | with kinetic 
energy, although the models with the highest mean absorption ve- 
locity all have ^kin > 1.45 x 10^^ erg. The lack of a general trend 
may seem contrary to expectations, but it reflects the impact of the 
different distributions of ignition points in each model series on the 
distribution of intermediate-mass elements in the ejecta. Within a 
given model series (i.e. for a given ignition setup), the mean |i;abs| 
increases with the kinetic energy of the explosion, which is di- 
rectly related to the criterion for deflagration-to-detonation transi- 
tion (^kin decreases steadily from dc5 to del within a model series; 
see Table [T]). We illustrate this for the DD2D_iso_06 model series 
in FiQ.\V7\( dotted line). 

In a recent paper, iFolev & KasenI (120 llh noted a correlation 
between the Si II 16355 absorption velocity at maximum light with 
^max _ y max pseudo-colouiQ at thc samc age, redder SNe la hav- 
ing higher |^;abs |. The lower panel of Fig. [18] shows the relation of 



^ difference between the B magnitude at B-band maximum and the V 
magnitude at V-band maximum. 



5.1 



mean rlap averaged over all viewing angles 

5.9 



6.7 



7.5 



8.3 



9.1 



9.9 






-16 
-15 
-14 
-13 
-12 
-11 



DD2D_asym_01_dc3 



DD2D iso 03 del- 



DD2D iso 08 dc3 



I - I f 

O 
■ DD2D_iso_06 model series (dc5 -► del) 



1.2 



1.3 



1.4 



1.5 



1.6 



Asymptotic kinetic energy [10 erg] 



Figure 17. Absorption velocity of the Si ii 16355 line at maximum light v^. 
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'^abs,max(Si II 16355) with intrinsic B — V colour at B-band max- 
imum (as opposed to B™^^ — ym^x pseudo-colour) in data from 
the Cf A SN Program. Although the correlation between both quan- 
tities is rather weak (r = —0.37; but note the large error on intrin- 
sic 5 — F c olour), the slope of th e relation is in agreement with 
that found bv lFolev & KasenI (1201 ih . For the 2D models of KRW09 
(Fig. [18] upper panel), the general trend is in the opposite direction 
(dashed black line). Within each model, the trend of absorption ve- 
locity with colour generally also has the opposite trend as that seen 

in th e data. 

iFolev & KasenI (1201 ll) provide a simple explanation for the 
observed correlation by invoking the greater line opacity in the 
5-band for SNe la with high- velocity ejecta, both through in- 
creased line overlap and a greater impact of Fe ll/Co II lines in 
the outer cooler regions of the ejecta where higher- velocity spec- 
tral features are formed . Based on radiative transfer calculations by 
iKasen & Plewal (120070 . they show tha t the detonating failed defla- 
gration (DFD) model of lPlewal (12007 ) is consistent with the asso- 
ciation of redder B — V colour with larger Si II 16355 blueshifts, 
when viewed from different viewing angles. As mentioned earlier, 
the 2D models of KRW09 also display variations of the Si II 16355 
absorption velocity with viewing angle, with associated variations 
in B — V colour, but the correlation in most cases is in the oppo- 
site direction as that se en in the DFD model, and in all cases much 
weaker than shown by iFolev & KasenI (l201ll) (the ir Fig. 8). Since 
the radiative transfer code used by KRW09 and Kasen & Plewal 
(2007) is the same, this discrepancy must be traced back to differ- 
ences in the explosion models (whether intrinsic or resulting from 
different physical approximations an d numerical treatment thereof) 
or nucleosynthetic post-processing jKasen & Plewal l2007l used an 
approximate alpha-network scheme which for instance does not in- 
clude sodium, while KRW09 used a more elaborate nuclear net- 
work to compute detailed abundances for a representative model, 
and interpolated the output for all other models). Analysis of more 
DFD models would be needed to check whether the correlation be- 
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Figure 18. Absorption velocity of the the Si ii 16355 line at maximum light 
v^". intrinsic B — V colour at S-band maximum, in the 2D models of 
KRW09 (upper panel) and in data from the CfA SN Program (lower panel). 
The colour-coding is the same as in Fig.|9] In both cases, we show a linear 
fit to the entire sample (dashed black line) and report the Pearson correla- 
tion coefficient (r) of ^abs,max with intrinsic B — V colour. We also show 
the linear relations obtained from fits to individual models (dashed colour 
lines). 



tween Si 1X16355 absorption velocity and maximum-light B — V 
colour is a generic feature of this type of explosion. This is ap- 
parently not the case for delayed- detonation models of SNe la, al- 
though Fig. [18] once again shows that the 2D models of KRW09 
lie in a different region of parameter space than the data (larger 
Si II 16355 blueshifts and redder B — V colours). 

One possible explanation resides in the different distributions 
of iron-group and intermediate-mass elements in both DFD and 
delayed-detonation models. In the DFD model, the deflagration is 
ignited on one side of the ejecta, but detonates on the opposite 
side, causing both IGE and IMF to be ejected at high velocities on 
the ignition side. In delayed-detonation models, the IGF and IMF 
are preferentially ejected in opposite directions (see Fig.[T] bottom 
row). Since the distribution of IGE sets t he B — V colour to a large 
extent (see, e.g. lKasen & Wooslevll2007h . one might expect the two 
detonation models to show opposite trends for |i;abs| and B — V. 



5.2 Spectral evolution 

The comparison of maximum-light spectral properties in the pre- 
vious section has several practical advantages: more data are avail- 
able at this age and they are of better quality (the use of a Monte 
Carlo radiative transfer code means this is also true of the model 
spectra, which consist of a higher number of photon packets at 
maximum light). Any thorough validation of a SN la model should, 
however, include a comparison of the time evolution of its spectra 



with observations. As the supernova expands, the layers where the 
spectrum is formed recede to deeper regions of the ejecta, where 
the composition is different and the expansion velocity smaller. By 
measuring various parameters of individual spectral features and 
their evolution with time, one has a complete census of a model's 
failures and successes: a poor model might reproduce certain spec- 
tral features at certain times, but a good model should reproduce all 
features at all times. 



5. 2. 1 Overall evolution 

In Fig. [19] we show the spectral evolution of model 
DD2D_iso_06_dc2 viewed along = 88° (see Fig. [3] for the 
corresponding UBVRIJHKs light curves), in five-day incre- 
ments between —10 and -1-20 d from B-band maximum (i.e. 
between ~ 10 d and ^-^ 40 d past explosion), compared to observa- 
tions of SN 2003du. Both models and observations correspond to 
aSNIawith Ami5(B) ?^ 1. 

At — lOd, the model spectrum displays a much bluer (hot- 
ter) SFD with weaker features than SN 2003du. The absorption 
due to Ca llH&K (^ 3700 A) is much weaker, and that due to 
the Ca II IR triplet (^ 8000 A) is non-existent, both clear signs of 
over-ionization in the ejecta. The weak Ca llH&K line certainly 
contributes to the excess of [/-band flux inferred from comparisons 
ofU — B colour curves (Fig.[T2l), but most of the discrepancy ap- 
pears to originate blueward of this, in a region dominated by ab- 
sorption by iron-group elements (grey highlighted region). This is 
most clearly seen in the -1-5 d spectrum of SN 2003du, which un- 
like the earlier ones extends blueward of 3600 A. As noted earlier 
based on the U — B colour curves, this excess of [/-band flux re- 
mains important up until -1-5 d, at which point the absorption due 
to Ca II H&K has significantly strengthened in the models though 
remains weaker than in SN 2003du. The S II doublet (-- 5400 A), 
well fit at —5 d is too weak in the models between -i-O d and -I-IO d. 
The same is true of the two iron-group dominated absorption com- 
plexes at ^ 4300 A and ^ 4800 A, which also reveal some of 
the model shortcomings from maximum light onwards. The ve- 
locity offset in the Si II 16355 absorption (broader yet shallower 
in the models; red dotted line) also becomes apparent at maximum 
light. Interestingly, the Si II 16355 absorption velocity appears to 
remain almost constant between — 10 d and -I-IO d while it decreases 
steadily in SN 2003du. 

From -1-10 d onwards, we highlight several discrepancies in 
other wavelength regions. The model appears to lack an absorption 
around -- 5000 A, which in SN 2003du and other SNe la is part 
of a complex Fe ll/Co II absorption feature. The pseudo equivalent 
width of this feature at +10 d is systematically ^-^ 50 A smaller in 
all the models with respect to observations. It is difficult to dis- 
entangle the effects of temperature (if too hot, as clearly the case 
before maximum light, this would delay the III^II recombination 
timescale of iron-group elements), the impact of the LTF approx- 
imation (non-LTF effects are particularly ini portant in treating the 
Fe ll/Co II opacities; see iBaron et al: 1996), and the high ejecta 
velocities (as clearly seen from the large blueshift of the Ca II IR 
triplet absorption)- which would enhance line overlap and smear 
out small-scale absorption features. 

Another striking feature is the lack of emission in the model 
around 5800 A. This feature is commonly attributed to the Na I D 
doublet. The absence of a strong 5800 A feature is a generic feature 
of these models, and not just particular to model DD2D_iso_06_dc2. 
The radiative transfer calculations we have at our disposal do not 
extend beyond -1-25 d or so past B-band maximum, so we cannot 
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Figure 19. Spectroscopic evolution of model DD2D_iso_06_dc2 viewed 
along 9 = 88° {red), compared to spectra o f SN 2003du (black; 
lAnupama. Sahu & Josell2005l:IStanishev et al.ll2007h . The synthetic and ob- 
served spectra have been scaled to match their respective rest-frame abso- 
lute B-band magnitude (see Fig. [3). The age of each spectrum is indicated 
in days from S-band maximum in the top-right corner of each plot. The 
grey shaded regions highlight the model tZ-band excess (between — lOd 
and -1-5 d), and several disagreements around 5000 A, 5800 A, and 8200 A 
(respectively attributed to Fe Ii/Co ii, Na I, and Ca ii) from -I-IO d onwards. 
The vertical red dotted line indicates the wavelength location of maximum 
absorption in Si ii 16355 in the models, which remains almost constant be- 
tween - 10 d and -1-10 d. 



say whether this line appears at later ages once the ejecta have 
cooled down, although KRW09 present synthetic spectra of model 
DD2D_iso_06_dc2 at +31 d that do not show evidence for this fea- 
ture (their Fig. 2). Artificially removing this line in the spectra of 
SN 2003du can lead to i?-band magnitude differences of up to 
~ 0.2 mag, and could explain part of the discrepancy in B — R 
and V — R colours between models and observations past maxi- 
mum light (see Fig.O. 

Non-LTE calculations by iBaron et alJ (120061) fail to repro- 
duce this feature (their Fig. 4), so the LTE approximation used by 
KRW09 is not necessarily the cause of this discrepancy. If this line 
is Na I D, the LTE approximation should be able to partly repro- 
duce this strong resonance line if sodium is abundant in its neu- 
tral stage (^Na = 10~^ throughout the ejecta in all the mod- 
els of KR W09, roughly corresp onding to one third of the solar 
value; e.g.. [Asplund et al.ll2009h . The absence of this line would 
then be related to an ionization problem (Na I has a low ioniza- 
tion potential of ^ 5.1 eV), which would likely persist until later 
times due to nont hermal ionization in the ever-thinning ejecta (e.g., 
iKozma et al.l 12005) . Another possibility is that this feature is not 
due to Na I, since its presence is difficult to reconcile with the pres- 
ence of singly- and doubly-ionized species, but we do not explore 
this possibility further at this stage. 



5.2.2 Absorption velocities 

Studying the blueshifts of individual absorptions and their evolu- 
tion with time enables one to indirectly probe the dynamics of the 
expansion. As a general rule, the blueshift decreases with time (as 
does |t'abs|), since the spectra form in deeper layers of the ejecta 
(i.e. at lower mass coordinate) as the SN expands. Various inho- 
mogeneities and radiative tra nsfer effects can, however, aff e ct this 
simple picture. For example, lOuimbv. Hoflich & Wheeled (l2007l) 
invoke a shell-like density structure (formed within the context 
of a pulsating delayed detonation or WD-WD merger scenario) 
to explain the ^^ 10 d-long plateau in absorption velocity of the 
Si II 16355 li ne in SN 2005hj. As for radiative transfer effects, 
IPessart et al.l (1201 1.) recently showed how taking into account time- 
dependent terms in the gas energy and statistical equilibrium equa- 
tions can lead to an increase of the He 1 110830 blueshift with time 
(contrary to all other He I lines) in the context of Type Ib/c super- 
novae. 

In what follows we explore the evolution of the absorption 
velocities in two strong lines in SN la spectra: Ca II 13945 (i.e. 
Ca II H&K) and the characteristic Si II 16355. The Ca II line profile 
spans a wavelength region where the doniinant source of opacity is 
due to bound-bound transitions (see, e .g. Pinto & Eastman 2000bl), 
and line overlap is an issue (see also iBlondin et al . 2006). More- 
over, observations of the Ca II IR triplet in early-time spectra of 
SNe la have revealed the pr esence of high- velocity absorption com- 
ponents fsee lMazzali et alii2005) . In what follows, and in the pres- 
ence of Ca II 13945 profiles with multiple absorptions, we report 
the velocity associated with the deepest absorption. The Si II 16355 
line is less affected by line overlap (at least until -i-lOd or so past 
B-band maximum), although some SNe la appear to have a high- 
velocity co mponent associated with this line at ea rly times (e.g., 
SN 2005cf: lGaravini et aLll2007l : IWang et aLll2009bh . 

In the top row of Fig. [20l we show the evolution of the 
Ca II 13945 absorption velocity, for four of our selected models 
(where we have again chosen model DD2D_iso_06_dc2 as repre- 
sentative of the DD2D_iso_06 model series), colour-coded accord- 
ing to viewing angle. The lower panel in each case shows simi- 
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lar measurements on SNe la in the same Amis (5) range from 
our spectroscopic sample, where the colour-coding is used to dis- 
tinguish individual supernovae. The models exhibit a strong di- 
versity in absorption velocity evolution. Model DD2D_iso_03_dcl 
has a steeply decreasing |i;abs| between — 15d and — 5d from B- 
band maximum (from 25000 km s~^ to 15000 km s~^), and lit- 
tle to no variation at later times. Such a variation is also appar- 
ent in some observed SNe la, although others display a steadier 
decrease in \vshs\ beyond -i-lOd. Model DD2D_iso_06_dcl {not 
shown', 0.85 ^ Ami5{B) ^ 1.08) has a similar behaviour. The 
other models of the DD2D_iso_06 series, however, all display little 
variation in the Ca II absorption velocity between — 15 d and nearly 
-1-10 d, with a sharp ^ 5000 km s~^ drop in I'^absl around -i-lOd, 
followed by a second phase of almost no variation. The discon- 
tinuous jump in Ca II absorption velocity could be in part related 
to a measurement artefact due to overlap and changes in relative 
strength between Si II 13858 and Ca II 13945. The same overall pat- 
tern is visible in model DD2D_iso_08_dc3, but it is not present in 
the data. Interestingly, model DD2D_asym_01_dc3 appears to dis- 
play both types of behaviour, depending on the viewing angle. 

As noted earlier, the Si II 16355 line is far less prone to mea- 
surement uncertainties, and is a better probe of the ejecta dynamics 
than Ca II 13945. The bottom row of Fig. |20l shows the evolution of 
the Si II 16355 absorption velocity in the same four models, again 
compared to observations in the same Ami^iB) range. Apart from 
model DD2D_iso_08_dc3, which displays little variation in |^;abs| 
between — 15d and maximum light, and a steady decrease there- 
after, the other three models show almost no variation for a given 
viewing angle (as expected, the largest scatter at any given age 
occurs in model DD2D_asym_01_dc3). Models DD2D_iso_03_dcl 
and DD2D_iso_06_dc2 even display a slight increase of |^;abs| with 
time. For most models we were unable to reliably measure the lo- 
cation of maximum absorption in Si II 16355 past -1-15 d or so. 

Also noticeable are the apparent lack of correlation between 
the Si II and the Ca II absorption velocities, as well as the ab- 
sence of high |^;abs| measurements before — lOd. Several SNe la 
from our spectroscopic sample display ^ 20000 kms~^ blueshifts 
at early times, where the models seem constrained to |^;abs| % 
15000 km s~^. Thus, while al l the model s are c onsistent with 
the "high- velocity" subclass of I Wang et al.l (12009a') based on the 
Si II 16355 absorption velocity at maximum light (see Fig.[T6]), they 
would be considered "normal" when considering the same mea- 
surement at — lOd. 



5.23 Physical interpretation and constraints on models 

Within the DD2D_iso_06 model series, model del shows a steady 
increase of |^;abs(Si II 16355) | with time, model dc2 shows almost 
no variation, and models dc3 through dc5 display a pre-maximum 
plateau followed by a progressively steeper post-maximum decline 
(Fig- ED right). The post-maximum velocity gradient is in the range 
50-150 km s~^ day~^ for these latter models, which represents a 
sizeable frac t ion of the range of velocity gradients observed by 
iBenetti et al.l (l2005h . The del through dc5 sequence corresponds 
to decreasing asymptotic kinetic energy within a model series (see 
Table [B, which directly affects the amount and radial distribution 
of intermediate-mass elements. 

The left panel of Fig. |2T] shows the angle- averaged silicon 
mass density profiles (= pXsi) for the dcl-dc5 models of the 
DD2D_iso_06 model series. Models with higher kinetic energy 
have a Si density distribution which peaks at higher velocities {r^ 
14000 km s"^ for model del cf. - 7000 km s"^ for model dc5). 



resulting in a higher \va,hs \ at any given age (see also Fig.[T7]). These 
models also synthesize more ^^Ni at the expense of intermediate- 
mass elements, hence both the height and width of the Si den- 
sity distribution decrease with increasing ^kin- This conditions the 
strength of the Si II 16355 line and its velocity evolution. In mod- 
els with higher kinetic energy, the lines are expected to be weaker, 
owing to the lower peak Si density, and to leave a more transient im- 
print on the spectra: in model del, the Si II 16355 absorption profile 
is identifiable out to ^ 15 d, while it is clearly visible out to ^ 20 d 
in model dc5. The narrower Si distributions (with less Si at low 
velocities) result in a more modest change of the absorption veloc- 
ity past maximum light: the Si II velocity gradient is almost flat in 
model dc2, and progressively steeper for models dc3 through dc5. 
The apparent increase in \vshs\ with time for model del is likely 
due to the weakness of the Si II line in that model, whose absorp- 
tion profile is biased by neighbouring lines. 

The absence of high- velocity (|vabs| ^ 15000 km s~^) mea- 
surements in the models at early times is a consequence of the 
drop in the silicon mass fraction in the outer regions of the ejecta 
{v > 15000 kms~^). To reproduce the Si II 16355 absorption pro- 
file in SN 2002bo, a SN la with a large abs orption veloc i ty at e arly 
times (^;abs ~ -18000 km s"^ at -13 d). IStehle et all (l2005h in- 
fer a homogeneous Xsi ~ 0.3 for v > 11000 km s~^. In model 
DD2D_iso_06_dc2, Xsi drops by two orders of magnitude (from 
- 0.3 to - 0.001) between - 16000 and - 21000 km s"\ The 
lack of a significant fraction of Si at high velocities in the models 
of KRW09 naturally accounts for the failure of the models to re- 
produce the high |^;abs| measurements observed in some SNe la, 
but we note that the nucleosynthetic yields in these low-density re- 
gions are subject to a large uncertainty. 

The absence of I'^absl > 25000 km s~^ measurements for 
Ca II 13945 is a consequence of the velocity cutoff used when 
remapping the hydrodynamical output on the 2D cylindrical grid 
for the radiative transfer calculations (see ^ 12.11 ). We do note, how- 
ever, that a low mass fraction Xca ~ 10 ~^ is sufficient to yield 
a Ca II absorption at the cutoff velocity, illustrating the strength of 
this doublet. 

The distribution of intermediate-mass elements typically ex- 
tends out to ^ 20000 km s~^ in the models of KRW09, regard- 
less of the initial distribution of ignition points. However, where 
almost all models with an isotropic distribution of ignition points 
have IMF down to the centre of the ejecta (Xime > 10~^ at 
V ~ Okms~^), most of the models with an anisotropic distri- 
bution (DD2D_asym) have no IMF below ^ 10000 km s~^, and 
hence display narrow 5000-10000 km s~^ "shells" of IMF in their 
ejecta. These IMF shells are a direct consequence of the weakness 
of the deflagration due to its one-sided ignition, resulting in a small 
pre-expansion and an IMF synthesis in a thin outer layer of the 
ejecta, and are clearly incompatible with the large > 10000 km s~ ^ 
widths of absorption profiles observed in SN la spectra. Not sur- 
prisingly, 5 of the 6 rejected models (see Table [B have narrow 
< 5000 km s"^ IMF shells in their ejecta. 



6 DISCUSSION: ARE ASYMMETRIC IGNITIONS 
RULED OUT BY OBSERVATIONS? 

In section 13.31 we rejected six models that do not yield satisfac- 
tory matches to observed SN la spectra for some viewing angles. 
Fig. I22I shows spectra for these models viewed along = 40° at 
maximum light, compared with observed maximum-light spectra 
in the same Ami^iB) range. The mismatch with observations is 
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Figure 20. Top row: Evolution of tlie Ca ii 13945 absorption velocity with time for four of our selected models {upper panels), colour-coded according to the 
viewing angle 6. The lower panels shows similar curves for SNe la in the same Ami5(B) range from our spectroscopic sample. The colour-coding in this 
case denotes individual SNe la, whose measurements are connected by dotted lines. Bottom row: same as above for the Si ii 16355 line. 
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Figure 21. Right: Angle-averaged silicon mass density (= pX^i) pro- 
files at f = 100 s past explosion for the DD2D_iso_06 model series (del 
through dc5). Left: IVIean absorption velocity curves for Si ii 16355 in the 
same model series. 



clear (and these objects would not be classified as SNe laQ but the 
synthetic spectra share some common properties: all have a very 
blue SED characteristic of hot ionized ejecta. Moreover, they show 
almost no sign of lines from intermediate-mass elements (Ca, S, 
Si) characteristic of SN la spectra, and instead are dominated by 



nor do they yield good matches to SN spectra of other types. 



iron-peak elements. This agrees fully with the output from the hy- 
drodynamical simulations: these models have the highest kinetic 
energy (^kin > 1.6 x 10^^ erg) and synthesize the most ^^Ni 
(M(^^Ni) > 1 M0), at the expense of intermediate-mass elements 
(M(IME) < 0.15 M0). The association of anisotropic distribu- 
tions of ignition points with greater amounts of synthesized ^^Ni 
holds for the other DD2D_asym models of KRW09. The defla- 
gration burns less material in these models, and by the time the 
detonation is triggered, the WD will have expanded less and the 
burning will occur at higher densities, increasing the fraction of 
material burnt to nuclear statistical equilibrium, and resulting in 
a globally symmetric ejecta. A direct observational consequence 
is that these explosions will tend to yield low Ami^iB) values. 
This is seen in this model series, where the DD2D_asym mod- 
els have Amis (5) values in a narrow range (Ami5(B) ^ 0.8), 
whereas the DD2D_iso models are evenly distributed across 0.9 < 
Ami5(B) < 1.4. 

These rejected models may cast some doubts on the valid- 
ity of the assumption of an anisotropic distribution of ignition 
points in SN la defla grations. This is in ap parent contradiction with 
the recent study by iMaeda et al.l (1201 Pal) , who associate the ob- 
served diversity in the spectroscopic evolution of normal SNe la 
with viewing angle effects in off-center delay ed-detonation models. 
By using nebular lines from stable iron-group ele ments as probe s 
of the distribution of deflagration ashes (see .Maeda et alj|2010bh . 
they showed that SNe la displaying low velocity gradients in the 
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Figure 22. Comparison of observed maximum-light spectra {black) with the 
subset of rejected models viewed along 6 = 40°, all in the range 0.76 ^ 
Ami5(B) ^ 0.89. SN1999gp and SN2005M are both 1991T-likeSNe la. 



Si 1X16355 line correspond to off-center explosions viewed from 
the ignition (deflagration) side, while SN with high velocity gra- 
dients resulted from the same explosion viewed from the opposite 
direction. 

In the present study, we see a variation of the Si II velocity gra- 
dient in explosions with an isotropic distribution of ignition points 
in which the criterion for deflagration-to-detonation transition is 
varied. This criterion primarily affects the asymptotic kinetic en- 
ergy of the explosion and results in a change in the distribution and 
abundance of silicon in the eie cta (see Fi^.[2T|). Wh ile this does not 
invalidate the interpretation of iMaeda et alj (l2010ah . it nonetheless 
shows that multiple parameters can affect the steepness of the Si II 
velocity gradient. 

We also note that the most DD2D_asym models of KRW09 
yield good matches to observations, such as DD2D_asym_01_dc3 
which is part of our subset of selected models. Moreover, not 
all DD2D_asym models with high explosion energies yield poor 
matches to observed spectra. Models DD2D_asym_07_dc2 and dc3 
have similar properties as the rejected models of Fig.[22l yet their 
maximum-light spectra are compatible with observed spectra of 
high-luminosity 1991T-like SNe la (see SN 1999gp and SN 2005M 



in Fig.ES, and with the peculiar SN 2000cx (Li et al.l l200ll) . both 
characterized by hot, ionized, and energetic ejecta. Asymmetric ig- 
nitions are thus not ruled out as a whole in the present study, but 
observations suggest that some anisotropic configurations do not 
occur in Nature. 

While the models with an asymmetric distribution of ^^Ni re- 
sult in a larger variation of photometric and spectroscopic prop- 
erties with viewing angle, we were unable to find specific mea- 
surements which correlate with the degree of asymmetry in the ex- 
plosion. Ejecta asymmetries are bou nd to leave an imprint o n the 
spectral line-profile morphology (see lPessart & HillieJuOllI in the 
context of Type II supernovae), but such signatures would probably 
be drowned in Monte Carlo noise in the synthetic spectra studied 
here. Spectropolarimetric observations of SNe la still provide the 
most direct (yet observationally exp ensive) probe of ejecta asym- 
metries (see IWang & Wheeled 120081 for a review), while nebular 
line profiles offer a means to assess the geometry of the explosion 
by probing the innermost regions of the ejecta. The limited number 
of photon packets used in the radiative transfer simulations and the 
assumption of LTE prevents us to investigate spectropolarimetric 
signatures or spectra during the nebular phase in the present study. 



7 CONCLUSIONS 

We have presented a detailed comparison of a recent survey of 
2D delayed-detonation explosion models by KRW09 with obser- 
vations of Type la supernovae. We apply standard methods used 
by SN la observers to compare the model light curves and spec- 
tra with empirical templates. This represents a significant step for- 
ward in the realism of the models. Running several light-curve fit- 
ters (MLCS2k2, SALT2, SNooPy) on synthetic {U)BVRI light 
curves, we find some tension between the light-curve shape of the 
models and actual data, the models having longer rise times. Based 
on cross-correlations with a library of SN la spectra, we quantified 
the overall resemblance of individual models to observed SNe la, 
and found that the best models/viewing angles lied systematically 
on the observed width-luminosity relation. 

Comparison of several photometric properties of the mod- 
els (rise times, maximum-light colours and their evolution with 
time) shows a broad agreement with observations, but reveals some 
problems with flux redistribution from the near-UV to the near- 
IR bands, a key mechanism needed to explain both the width- 
luminosity relation and the secondary maxima in the NIR light 
curves, and mediated by the lll ^ll reco mbination timescale of 
iron-group elements ( Kasen 2006 : iKasen & Woosle v 2007). Sub- 
sequent investigation of spectra for a subset of selected models 
confirmed the excess of [/-band flux in the models at early times, 
likely caused by a hot ionized ejecta and subsequent lack of ab- 
sorption by Fe ll/Co II. Interestingly, one of our selected models 
(DD2D_iso_08_dc3) that shows the best overall agreement in opti- 
cal and NIR colour evolution with observations lies off the width- 
luminosity relation (i.e. the colours match, but not the luminos- 
ity). This reveals one limitation of our approach, which relies on 
maximum-light spectra to rank the different models, whereas a 
combination of photometric and spectroscopic properties is needed 
for a proper evaluation. 

Comparison of maximum-light spectra show the models have 
systematically large absorption velocities (most visible in the 
Si II 16355 line), affecting the relative shapes and strengths of spec- 
tral features and smoothing out small substructures observed in 
iron-dominated absorption complexes at ^ 4300 A and ^ 4800 A. 
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Consequently, correlations between several spectroscopic indica- 
tors and Ami5{B) decline rate have a much larger scatter in 
the models. The relation found bv lFolev & KasenI (120111) between 
Si II 16355 absorption velocity and intrinsic B — V colour (redder 
SNe la having larger |i;abs|) is also not reproduced in the mod- 
els, most showing a weak correlation in the opposite direction. 
However, the observed correlation is weak and subject to a large 
uncertainty given the errors on intrinsic B — V colour inferred 
from the data. Nonetheless, we identify a trend of larger absorp- 
tion blueshifts for higher kinetic energy for models in which only 
the criterion for deflagration-to-detonation transition is varied. 

The overall evolution of the model spectra compares well with 
observations, as illustrated by the comparison of synthetic spectra 
for model DD2D_iso_06_dc2 between -lOd and -h20d from B- 
band maximum with observed spectra of SN 2003du, but several 
discrepancies characteristic of most models are apparent. The syn- 
thetic spectra are too blue (i.e. too hot) at early times, and the ion- 
ization that ensues affects their subsequent evolution. Most notably, 
the models appear to lack an absorption feature around ^ 5000 A 
(attributed to Fe ll/Co ll), and fail to reproduce the strong emission 
feature at ^ 5800 A (attributed to Na iD) from -hlOd onwards. 
Non -LTE calc ulations of SN la spectra also fail to reproduce this 
line tearon et al. 2006). 

The evolution of the Ca II 13945 absorption velocity with 
time exhibits a strong diversity in the models which contrasts 
with the steady and smooth decrease seen in the data. Model 
DD2D_asym_01_dc3 has an asymmetric ^^Ni distribution and 
shows all types of behaviour depending on the viewing angle, due 
to the varying radial Ca distribution with different lines of sight. 
The evolution of the Si II 16355 velocity shows little variation be- 
fore maximum light, while the post-maximum evolution (the veloc- 
ity gradient) appears conditioned by the kinetic energy of the ex- 
plosion, affecting the abundance and radial distributi on of silicon 
While t his does not contradict the recent findings of iMaeda et al.l 
(1201 Par) , who associate the observed diversity in velocity gradients 
with viewing angle effects in off-center explosions, it shows that 
the interpretation of such gradients depends on more than a single 
parameter of the explosion. 

We reject six models of KRW09 with highly asymmetric igni- 
tion conditions and are characterized by large amounts (> 1 M©) 
of ^^Ni. We do not reject off-center delayed-detonation mod- 
els for SNe la as a whole, but note the extreme sensitivity of 
the amount and distribution of burning products in the deflagra- 
tion phase to the initial dist ribution of ignition points (see also 
iLivne. Asida & Hoflichll2005l) . 

Throughout this paper we have focused on discrepancies be- 
tween the models and observations more than we have highlighted 
their mutual agreement, but this merely results from the unprece- 
dented level of detail of our study. Such detail is necessary to use 
the predictive power of the models to provide a physical basis to 
some observed trends, as well as use the data to impose meaningful 
constraints on the models. The 2D delayed-detonation models of 
KRW09 have a degree of fidelity which makes them amenable to 
the same analysis we use on observations of SNe la, but they still 
require some adjustments to accurately match the data. The ability 
to reproduce the bolometric/multi-band light curves and the width- 
luminosity relation is a necessary but not a sufficient condition for 
a model to be considered a valid approximation of real SNe la. Fur- 
ther insights from three-dimensional hydrodynamical simulations, 
more accurate nucleosynthetic post-processing, and full non-LTE 
radiative transfer calculations are all part of the solution. We are 
confident that a detailed comparison of light curves and spectra 



from grids of models using the framework developed in this paper 
will lead to a better understanding of SN la explosion mechanisms. 
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Table 1. Asymptotic kinetic energies, abundances, peak bolometric luminosities and decline-rate ranges for the 44 2D delay ed-detonation models of KRW09. 



Model 


^kin 


M(56Ni) 


M(stable IGE) 


M(IME) 


M(0) 


M(C) 


-^bol,peak 


Ami5(B) 


Notes 




(10^1 erg) 


(M0) 


(M0) 


(M0) 


(M0) 


(M0) 


(10^3 ergs-1) 






DD2D_iso_01_dc2 


1.548 


0.942 


0.199 


0.215 


0.044 


0.007 


1.921-2.133 


0.91-1.03 




DD2D_iso_01_dc3 


1.506 


0.889 


0.194 


0.253 


0.061 


0.011 


1.818-2.026 


0.99-1.08 




DD2D_iso_01_dc4 


1.463 


0.832 


0.193 


0.289 


0.078 


0.016 






not included 


DD2D_iso_01_dc5 


1.442 


0.809 


0.192 


0.302 


0.085 


0.019 


1.614-1.812 


1.07-1.17 




DD2D_iso_02_dc2 


1.569 


0.901 


0.256 


0.210 


0.034 


0.007 


1.863-2.172 


0.88-1.16 




DD2D_iso_02_dc3 


1.515 


0.796 


0.253 


0.298 


0.051 


0.010 


1.563-1.915 


0.97-1.25 




DD2D_iso_02_dc5 


1.410 


0.612 


0.241 


0.456 


0.082 


0.016 


1.281-1.430 


1.15-1.48 




DD2D_iso_03_dcl 


1.509 


0.799 


0.274 


0.280 


0.045 


0.009 


1.763-1.850 


0.79-0.97 


subset 


DD2D_iso_03_dc2 


1.385 


0.584 


0.263 


0.448 


0.097 


0.016 


1.277-1.357 


0.90-1.16 




DD2D_iso_03_dc3 


1.288 


0.441 


0.257 


0.548 


0.138 


0.028 


0.925-1.017 


1.08-1.33 




DD2D_iso_04_dcl 


1.502 


0.774 


0.275 


0.298 


0.051 


0.011 


1.499-1.938 


0.77-0.99 




DD2D_iso_04_dc2 


1.379 


0.562 


0.264 


0.470 


0.095 


0.017 


0.983-1.519 


0.85-1.19 




DD2D_iso_04_dc3 


1.263 


0.399 


0.255 


0.583 


0.145 


0.027 


0.633-1.228 


1.06-1.38 




DD2D_iso_04_dc4 


1.226 


0.369 


0.252 


0.593 


0.162 


0.033 


0.581-1.147 


0.99-1.44 




DD2D_iso_05_dcl 


1.470 


0.718 


0.260 


0.357 


0.060 


0.011 


1.560-1.594 


0.87-1.11 




DD2D_iso_05_dc2 


1.333 


0.496 


0.249 


0.530 


0.114 


0.019 


1.111-1.131 


1.06-1.22 




DD2D_iso_05_dc3 


1.210 


0.330 


0.237 


0.648 


0.160 


0.033 


0.754-0.826 


1.21-1.44 




DD2D_iso_05_dc4 


1.151 


0.293 


0.232 


0.648 


0.189 


0.046 


0.691-0.725 


1.23-1.44 




DD2D_iso_06_dcl 


1.552 


0.882 


0.256 


0.227 


0.036 


0.009 


1.859-1.959 


0.85-1.08 


subset 


DD2D_iso_06_dc2 


1.446 


0.698 


0.247 


0.378 


0.073 


0.012 


1.443-1.564 


1.06-1.32 


subset 


DD2D_iso_06_dc3 


1.371 


0.567 


0.240 


0.485 


0.099 


0.019 


1.151-1.288 


1.21-1.48 


subset 


DD2D_iso_06_dc4 


1.312 


0.472 


0.235 


0.561 


0.117 


0.025 


0.972-1.064 


1.30-1.53 


subset 


DD2D_iso_06_dc5 


1.275 


0.421 


0.231 


0.599 


0.129 


0.029 


0.896-0.959 


1.34-1.57 


subset 


DD2D_iso_07_dc2 


1.282 


0.444 


0.246 


0.556 


0.139 


0.023 


0.867-1.079 


1.04-1.24 




DD2D_iso_07_dc3 


1.181 


0.346 


0.238 


0.603 


0.183 


0.039 


0.640-0.861 


1.18-1.30 




DD2D_iso_08_dcl 


1.477 


0.719 


0.275 


0.341 


0.060 


0.012 


1.534-1.767 


0.80-1.07 




DD2D_iso_08_dc2 


1.302 


0.447 


0.259 


0.549 


0.130 


0.023 


0.930-1.072 


1.00-1.38 




DD2D_iso_08_dc3 


1.204 


0.329 


0.250 


0.624 


0.169 


0.035 


0.695-0.785 


1.16-1.39 


subset 


DD2D_asym_01_dc2 


1.416 


0.677 


0.224 


0.410 


0.081 


0.015 


1.095-1.811 


0.93-1.09 




DD2D_asym_01_dc3 


1.379 


0.644 


0.221 


0.419 


0.103 


0.021 


0.997-1.696 


0.96-1.15 


subset 


DD2D_asym_02_dc2 


1.387 


0.614 


0.229 


0.454 


0.096 


0.016 


1.234-1.451 


0.96-1.23 




DD2D_asym_02_dc3 


1.252 


0.458 


0.216 


0.558 


0.147 


0.027 


0.886-1.139 


1.18-1.38 




DD2D_asym_03_dc2 


1.580 


0.982 


0.234 


0.157 


0.028 


0.006 


1.895-2.460 


0.79-0.88 


rejected 


DD2D_asym_03_dc3 


1.572 


0.965 


0.233 


0.171 


0.031 


0.007 


1.834-2.426 


0.79-0.88 


rejected 


DD2D_asym_04_dc2 


1.577 


0.974 


0.201 


0.190 


0.036 


0.007 


1.722-2.262 


0.80-0.91 




DD2D_asym_04_dc3 


1.562 


0.955 


0.198 


0.204 


0.042 


0.009 


1.666-2.243 


0.81-0.96 




DD2D_asym_05_dc2 


1.538 


0.939 


0.208 


0.211 


0.042 


0.007 


1.979-2.054 


0.77-1.00 




DD2D_asym_05_dc3 


1.508 


0.883 


0.205 


0.253 


0.057 


0.011 


1.876-1.981 


0.82-1.06 




DD2D_asym_06_dc2 


1.628 


1.082 


0.221 


0.087 


0.014 


0.003 


2.175-2.562 


0.78-0.88 


rejected 


DD2D_asym_06_dc3 


1.622 


1.068 


0.213 


0.104 


0.019 


0.004 


2.144-2.537 


0.78-0.89 


rejected 


DD2D_asym_07_dc2 


1.579 


1.053 


0.175 


0.144 


0.030 


0.005 


2.192-2.362 


0.74-0.94 




DD2D_asym_07_dc3 


1.558 


1.034 


0.170 


0.158 


0.038 


0.008 


2.164-2.332 


0.76-0.95 




DD2D_asym_08_dc2 


1.632 


1.103 


0.206 


0.081 


0.014 


0.003 


2.113-2.586 


0.76-0.87 


rejected 


DD2D_asym_08_dc3 


1.617 


1.081 


0.189 


0.113 


0.020 


0.005 


2.043-2.578 


0.77-0.89 


rejected 
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